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Abstract 

We present a general formalism able to derive the kinetic equations of polymer 
dynamics. It is based on the application of nonequilibrium thermodynamics 
to analyze the irreversible processes taking place in the conformational space 
of the macromolecules. The Smoluchowski equation results from the anal- 
ysis of the underlying diffusion process in that space within the scheme of 
nonequilibrium thermodynamics. We apply the method to different situa- 
tions, concerning flexible, semiflexible and rod-like polymers and to the case 
of more concentrated solutions in which interactions become important. 

I. INTRODUCTION 

In the theory of polymer dynamics ||, 0, and in general in the dynamics of complex 
fluids, the knowledge of the kinetic equations describing the evolution of the probability 
density for the configurations of the suspended phase plays a central role. The formulation 
of Fokker-Planck or Smoluchowski equations constitutes the starting point in the study of 
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the configurations of the system and in the characterization of its global behaviour. It is then 
of primary importance to establish methods able to provide those equations for a complex 
system possessing a huge number of interacting degrees of freedom and exhibiting a wide 
variety of time and length scales. 

^From the statistical mechanical point of view, a polymer solution can be conceived as 
a set of objects of mesoscopic size immersed in a heat bath. The different configurations 
adopted by the macromolecules depend on their intrinsic nature, the interactions between 
their different parts and on the interactions with the bath. All these ingredients must be 
considered in the formulation of Fokker-Planck equations. 

Our purpose in this paper is to present a general and simple formalism able to derive 
kinetic equations of the Fokker-Planck type for polymers in solution. The formalism is based 
on the fact that processes leading to variations in the conformation of the macromolecules 
can be described by nonequilibrium thermodynamics. The extension of that theory to the 
mesoscopic level of description, referred to as mesoscopic nonequilibrium thermodynamics 
(MNET), has been applied to different situations pertaining to the domain of transport 



phenomena ||- || and activated processes ||, [TCfl . We will show that it can also be used 
in the case of polymers to systematically derive the Smoluchowski equations for polymers 
of different nature. 

The paper is organized in the following way. In section 2, we present the general deriva- 
tion of the kinetic equation for a macromolecule in solution in the framework of MNET. 
These results are applied in Section 3 to the case of flexible and rod-like polymers to derive 
the corresponding Smoluchowski equations. The case of semiflexible polymers is presented 
in section 4. We derive the functional Fokker-Planck equation for the probability of the 
different configurations of the chain. In Section 5, we treat the case of more concentrated 
solutions in which interactions between segments belonging to different chains become im- 
portant. The functional Fokker-Planck equation is obtained in terms of the concentration 
of segments. Finally, in Section 6 we present our main conclusions. 
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II. KINETIC EQUATIONS FOR THE DYNAMICS OF MACROMOLECULES IN 

SOLUTION 



A polymer solution is usually modeled as a suspension of macromolecules in a Newtonian 
solvent [|TJ , [TT| . This simplification may be carried out due to the existence of well-separated 
time and length scales related to the carrier fluid and to the macromolecules. The dynamics 
of the fluid itself is then governed, at the continuum level, by the Navier-Stokes equations 
whereas the polymers are considered as suspended objects whose presence has implications 
in the dynamics of the whole system. Since the interaction between the polymer and the 
solvent occurs through its contour, the form of the polymers will be of crucial importance 
in the characterization of the properties of the system. 

The polymer solution is usually conceived as a two-component system. One of them, is 
the host fluid that plays the role of a heat bath at constant temperature T. To characterize 
the different configurations of the macromolecule, we will introduce the vector 7 whose 
components constitute the independent parameters necessary to completely determine its 
configuration. We can then define the quantity P( / y,t)/N as the probability density for 
the polymer in the configuration space (7-space), where N is the number of polymers in 
suspension. 

Since P(j, t) is a conserved quantity, it evolves according to the continuity equation 

— = -V 2 - J(7,t) (1) 

with J(7, t) being a current in the 7-space. 

Non-equilibrium thermodynamics establishes the Gibbs equation in which changes in the 
entropy are related to variations in the density of the particles. We will generalize that idea 
by assuming that variations in the probability density in 7-space are responsible for changes 
in the total entropy of the system, S. Our starting point will then be to formulate the Gibbs 
equation 

SS = ~J //(7,t)5P(7,^7, (2) 
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where 7/(7, t) is a chemical potential defined in 7-space. This equation holds for constant 
energy and volume of the system. 

The expression of the chemical potential can be identified by imposing that the Gibbs 
equation (0) must be compatible with the Gibbs entropy postulate. This postulate estab- 
lishes that 

S = -k J P In P/P Leq -d 1 + S Leq - , (3) 

where k is the Boltzmann constant, and the local equilibrium distribution p Le i- is 

P^-(7, t) = exp{(iV yUp - U) /kT} , (4) 

with fi p being the chemical potential of the polymers at local equilibrium, and U the potential 
energy. In addition, S q ~ is the total entropy at local equilibrium whose variations are given 
by the Gibbs equation 

8S 1 ^ = ~ix p 5N . (5) 

By taking variations in Eq. (|]) and comparing with Eq. (Q), after using Eq. (|5]) we obtain 

the expression for the chemical potential /x(7, t) 

kT kT U 

M7, t) = n P + — ln(P/P^-) = — In P + - . (6) 

The Gibbs equation then resembles its corresponding expression for an ideal mixture in 
which the different species would correspond to the different values of 7, which may be 
interpreted as an internal coordinate or degree of freedom |J. 

The entropy production, accounting for dissipation in the 7-space related to variations 
of the configuration of the polymer follows straightforwardly from the Gibbs equation (0). 
Taking the time derivative in Eq. (^|) and using Eq. (|I|), after partial integration we obtain 

<r=~ /j(7,t)"V 2 M7,^7- (7) 

Following the rules of nonequilibrium thermodynamics || , from the entropy production one 
derives the linear laws 



J(7, t) = ~ J Hi, i f )VMi f , t)dij , (8) 

where L(j, 7/) is a matrix of phenomenological coefficients. Its local version, occurring when 
L(7, 7/) = L(j)5(j — 7/), is then given by 

Z(7,t) = -^(7)V 2 /i(7^)- (9) 

Inserting the expression of the chemical potential Eq. @ into Eq. (|9|) and the resulting 
current into the continuity equation (JJ) one obtains 

or 

— = V 2 • M ■ (kTV x P + PV 2 U) , (II) 

where use has been made of the relation between the matrices of phenomenological coeffi- 
cients and the matrices M 

M=-^. (12) 
— NTP v J 



If we interpret M as the mobility matrix related to the dynamics of the vector 7, Eq. (|Tl~D 
is precisely the Smoluchowski equation describing the evolution of the probability density 
in 7-space. To solve the Smoluchowski equation, we need to know the expression for the 
mobilities. This quantities can be borrowed from hydrodynamics |Oj, [[13 . 



III. FLEXIBLE AND ROD-LIKE POLYMERS 

As a first example, we will apply the results of the previous section to study the dynamics 
of a dilute solution of flexible polymers. The macromolecules can be modeled as a set of 
n segments (monomers) whose positions will be denoted by 7 = (Ri, R 2 , R n ) PJ- The 
quantity 7 is then considered as the internal coordinate and its different values are related 
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to a point in a 3n- dimensional configurational space. We can now introduce the density 
P(7, t) which evolves according to the continuity equation 

f)P n 

ir-Ev r (Ji + Pvi) , (13) 
UL i=l 

where Jj is the current associated with the diffusion of the probability density along the 
coordinate Ri, and Vi(Ri,t) = vo(r = Ri,t), with vq being an externally imposed velocity 
field. 

We can also define the Gibbs equation from Eq. (^) for the polymer phase. As we did in 
the previous section, from that equation, and using Eq. (fl3D it is possible to compute the 
entropy production related to the diffusion process in the conformational space 

^ = -i/E^'V A A*d{fli}. (14) 
1 J i=i 

From this quantity one then derives the linear laws 

-* 1 n 

Ji = -TpY.^ij-^Rit 1 ' ( 15 ) 
1 3=1 

where are now matrices of phenomenological coefficients which according to the Onsager's 
reciprocity principle satisfy the relations 

hi = 4 ■ ( 16 ) 

Here the superscript t stands for the Hermitian adjoint. The chemical potential is also given 
by Eq. (H), where now the potential energy U{^,t) may in general include contributions 
due to intermolecular interactions among monomers, external forces or even contributions 
coming from the presence of obstacles. The Smoluchowski equation then follows from Eqs. 
(0), (0), and ©: 

BP n 

¥ = EV r Mij ■ (kTV M P + PV Mj U) -£V 4 - (Pi*) • (17) 

As a second example we will consider the case of a rod-like polymer ]]J which will be 
modeled as a cylinder of length L and diameter b. The internal coordinate is now taken as 
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7 = (R, ft), with R being the centre of mass position vector and n the unit vector pointing 
along the axis of the polymer. 

The continuity equation for P(R, ft, t) is now written as 

dP - 

= ~Vr ■ Jr - ■ Jn , (18) 

where and Ja are the corresponding diffusion currents. Following the steps indicated 
before we then obtain the entropy production which has the usual form of a sum of fluxes- 
force pairs 

a = — — J dRdn ■ V^u — — J dRdn Ja ■ Vafi . (19) 
From this equation one then derives the linear laws 

(20) 













Jr = 


rp^-R,R 






> n ■ V ^ 
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V^/i - 


rp—n,. 





where L^ B , a, (3 = (R, ft), are matrices of phenomenological coefficients, written in local form. 
Alternatively, by using the expression for the chemical potential, /z(n, R,t) = ^ In P(n, R, t) 
(for simplicity sake we do not consider the action of any field of forces) , these equations can 
be rewritten as 

Jr = -Ma* ■ VrP - M R>a ■ V H P (21) 

where use has been made of the relation among the mobilities and the phenomenological 
coefficients. According to Onsager's symmetry principle, the phenomenological coefficients 
satisfy the relations 

Ma- = Aft 3 . (22) 

R,n fi t R v I 

By substituting the linear equations (^Tj) into the continuity equation ( PD we then obtain 
the kinetic equation 



Up to lowest order in the parameter b/L, it is known [JTJ that the mobility matrices can 
be written as 

M.RR = lA™ + A*l(l — nn) 
Kn,n = /i r (I - nn) 

% = M.Sfl = , (24) 

where fit ± and /i r are the translational and rotational scalar mobilities, and 1 is the unit 
tensor. Insertion of Eqs. fl24j) into (|23| ) then leads to the Smoluchowski equation 



OP 

— = V n ■ { [D\\nn + D ± (l - fin)] • V^P} + .D r 7? 2 P , (25) 

where we have introduced the diffusion coefficients D\\ = kTfA, D± = kTf/ ± , and D r = kTfi r 
and TZ = n x d/dn is the rotational operator. 

IV. SEMIFLEXIBLE POLYMERS 

Unlike flexible polymers, semiflexible polymers characterize by the presence of some 
degree of stiffness that introduces a persistence length []1~4| |. The semiflexible polymer may 



be described by a continuous curve r(l), with I (0 < I < L) being the contour length along 
the chain. To derive the Fokker-Planck equation, we will introduce the continuity equation 

^n{m,t) = -jdl'^ryJ{l',t). (26) 

Here A r_1 P({r(/)}, t) represents the probability density for the configuration f(l) at time t, 
J{l,t) is the diffusion current in the configuration space. 
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The rate of change of entropy can be obtained from the corresponding Gibbs equation, 
and reads 



J 5f(lM{f(l)ht)^Pmi)ht) ■ (27) 



dS 1 r .s d 

~dt ~~f 



Here, the chemical potential is given by 

kT 1 
MM/)},*) = ^lnP({r(/)},t) + ^£/({r(0» , (28) 

where the potential U({r(l)}) consists of entropic and elastic energy of bending contributions 



«'>4-r(f)-i/:(f)^' 

with Iq being the effective segment length, and e the bending force constant. 

By making use of Eq. fl26|), and after partial integration we obtain the entropy production 
related to the diffusion process in the configurational space 

,= jdlJ-^. (30) 

From the entropy production we derive the linear laws 



J(l,t) = ~ J dl' L(l,l') ■ , (31) 



where we have taken into account that now locality is not fulfilled. The matrices of Onsager 
coefficients L(l, I') satisfy de Onsager's relations 

L(M') =L{l',iy . (32) 

Using Eqs. (|31| ) with ( ffip in the continuity equation (|26|), we obtain the functional Fokker- 
Planck equation 

f = / dl I dl W) ■ * l < f ' ■ { kT W) p + p m um ' )]) } (33) 

describing the evolution of the probability density for the configurations of the polymer. 



V. MANY CHAIN SYSTEMS 



In this section, we will indicate how the formalism developed in previous sections for 
dilute polymer solutions also applies to the case of higher concentrations. We will analyze 
the case of flexible polymers and consider as the internal coordinate the vector of the con- 
figuration space 7 = {R„}, with i?" being the position vector of the v-th segment pertaining 
to the a-th chain. The probability density in the configuration space is then N~ 1 P({R^}, t) 
where now N must be understood as the total number of segments in suspension. As we did 
in section 3, the starting point in our analysis is the partial Gibbs equation for the polymer 
phase. This equation together with the continuity equation 

f = -EV^J: (34) 

allows us to express the entropy production as 

ff = -^/E^-V^4K}, (35) 

where J" are currents in the configurational space spanned by the values of the coordinates 
{R„}- Moreover, the chemical potential consists also of ideal and non-ideal contributions 
and is given by 

kT -, 1 

KWh t ) = Jf l * t) + j^U({Rl}) . (36) 

The linear laws are then written as 

^ = -^E^"V^, (37) 

b,r] 

where are matrices of phenomenological coefficients relating the current J", associated 
with the internal coordinate {R^,}, with its conjugated thermodynamic force V^/i. These 
matrices satisfies the Onsager relations 

rab rfeat foo\ 

LLvrf — Ll V v ■ K 6 °) 

Substitution of Eq. ([37|) in Eq. (|34] ) then leads to the kinetic equation 
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BP 

-qT = ~ E • M-tr, ■ {kTV nb P + FVj^U} . (39) 

a, b.u r] 

Alternatively, we could now proceed to describe the distribution of the polymer segments 
in terms of the local segment density c(r) = J2a,v$(r — R%) This description, which is 
useful when performing ensemble averages of quantities depending on the local segment 
density, would then imply the introduction of the probability density P({c(k)}, t) ([HJ in the 
Gibbs equation in order to specify the local state of the polymer solution. Proceeding along 
the lines indicated in previous cases, we obtain the entropy production 

° = -\\ 5c{k) J dkJ({c(k)},t;k)^^{c(k)},t) , (40) 

where J({c(k)},t; k) is the corresponding current and the chemical potential is given by 

kT - 1 

A*({c(*)}) = -jy M{c(*)}) + j^U({c(k)} . (41) 

In Eq. (^UD we have introduced functional integrals and derivatives because the continuum 
character of the index k, with k z > 0. The linear relations are now written in non-local form 

as 

J({c(k)},t;k) = -U dk'L(k,k')-^u({c(k'),t) . (42) 
1 J dc(k') 

Following the procedure indicated previously we finally arrive at the Fokker-Planck equation 



^-P({c(k)}) = [die f dk'—^M(k, k') \ kT—^P({c(k')\) + P-L-U({c(k')}) 
dt J J Sc(k) I 6c{k') Sc(k') 



(43) 



where for the sake of simplicity we have not considered the effects of an external flow. 
For the case in which the matrix L(k, k') is local in Fourier space due to the preaveraging 
approximation, this equation agrees with the one obtained in Ref. |l] from the Langevin 
equation for the density. In the method we have presented, the kinetic equation comes from 
the analysis of the dissipation in phase space, according to the principles of nonequilibrium 
thermodynamics. 
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VI. CONCLUSIONS 



In this paper, we have applied the method of mesoscopic nonequilibrium thermodynamics 
to derive the kinetic equations governing the dynamics of macromolecules in solution. The 
Fokker-Planck equations can be obtained from the analysis of the underlying driven diffusion 
process occurring in the space of configurations of the polymers. 

The intrinsic nature of the polymer defines the set of degrees of freedom necessary for its 
complete dynamic description. One then formulates the continuity equation for the probabil- 
ity density in the configuration space which introduces an unspecified current in that space. 
The determination of that current is carried out through nonequilibrium thermodynamics 
giving rise to the Fokker-Planck equation. The simplicity of the procedure contrasts with 
the more elaborated method of deriving the Fokker-Planck equation from the correspond- 
ing Langevin equation, becoming more cumbersome for complex systems involving many 
degrees of freedom. 

The formalism we have proposed can be applied to polymers of different nature in the 
presence of external gradients |16[ or forces and for arbitrary values of the concentration. It 
may then constitute a useful tool in the theory of polymer dynamics. 
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